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> ! Abstract 
CO 

The Hamiltonian structure of the two-dimensional dispersionless Toda hierarchy is stud- 
' ied, this being a particular example of a system of hydrodynamic type. The polynomial 

conservation laws for the system turn out, after a change of variable, to be associated 
with the axially symmetric solutions of the 3-dimensional Laplace equation and this 
enables a generating function for the Hamiltonian densities to be derived in closed form. 



o 



1. Introduction 



In this paper the Hamiltonian structure of the integrable system 

S(x,t) t = P(x,t) x , 
P(x,t) t = P(x,t)S(x,t) x , 

(or equivalently, the equation (\ogP) tt = P xx ) will be investigated. The reason for the 
interest in this system is two-fold. Firstly, it appears in various physical applications, for 
example: 

• In the construction of self-dual Vacuum and Einstein- Weyl metrics [0. 

• In topological field theory, as a solution to the WDVV equations 0. 

• As the long- wave (or dispersionless) limit of the Toda equation |3|, [|. 

Secondly, despite its simplicity, it has the same generic properties as more complicated 
systems, so it can be used as a 'test bed' for various ideas. (See for example the recent 
paper of Mineev-Weinstein[^|.) It is one of the simplest examples of an equation of 
hydrodynamic type that is, an equation of the form 

4 = v?(u)t4 (2) 
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about which there is an extensive literature || M. Much of this, however, involves general 
theory, and there are comparatively few specific examples. The purpose of this paper is to 
derive concrete results for the dispersionless Toda system. Rather than study ([[]) directly, 
a transformation will be made to a new set {u,v} of variables in which the equations 
are symmetric under the interchange of u and v , this enabling results from the theory 
of homogeneous symmetric functions to be used. Explicitly, let {S, P} be a solution to 
([!]) and let {u,v} be the roots of the quadratic 

z 2 - Sz + P = . 

In these new variables the original system ([!]) becomes 

u t = uv x , , . 

v t = vu x , 

Note that if u = v then this reduces to the dispersionless Korteweg-deVries equation. 

This change of variables may be used to see the following property of any solution to 
(0). Given {u,v} satisfying ([D a new solution is given by the roots of the quadratic 

z 2 — (u + v + k)z + uv = . 

Moreover, one has a version of the 'theorem of permutability' -applying this transfor- 
mation first with constant k% then with &2 results in the same solution as applying the 
transformation with constant &2 followed by that with k\ . The general solution to ([!]) or 
(^) may be found implicitly by performing a hodographic transformation. Specific exam- 
ples have been calculated in || and the general form of the solution in ||. 

Another motivation comes from the construction of integrable, relativistic system in 
(1 + 1) dimensions. Any solution to (Q) of the form 

f (a±a±\ 

J \dt dx) 

9_mi 

(or, equivalently, a similarity solution) will result in such a system for <p , with x and t 
being light-cone variables. The equations can be reduce to the single first order differential 
equation 

dY _ Y(l + 2X - XY) 
dX ~ X(1 + 2Y- XY) ' 

this being an Abel equation of the second kind. Despite the integrability of the original 
system (J^) we have not been able to solve this equation, except in the case Y = X (or 
u = v) which results in the well-known Born-Infeld equation for <fi 0. 

The rest of this paper is arranged as follows. In section 2 conservation laws will be 
derived for this system, together with a generating function for them, and in section 3 



u = 



v = 



the associated hierarchy will be studied. In section 4 these results will be used to study 
the Hamiltonian structure of this system and its hierarchy. These results follow from 
direct calculation and do not use the corresponding Lax pair. This Lax pair has a very 
simple solution, details of which may be found in the appendix. Various generalisations 
are presented in section 5. 



2. Conservation Laws 



By direct calculation it is easy to find the first few polynomial conservation laws: 



[u + v] t = [uv] x , 
[u 2 + Auv + v 2 ] t = [2uv(u + v)] x , 
[u 3 + 9uv(u + v) +v 3 ] t = [3uv(u 2 + 3uv + v 2 )] x , 

and a general such law will be denoted by 

Qn,t = Fn,x , (4) 

where Q n (the charge density) is an n th -order symmetric polynomial in u and v , and T n 
(the flux density) is an (n + l) th -order symmetric polynomial in u and v . The coefficients 
in Q n are the squares of the binomial coefficients B(n, r): 

n 

2„.r„,n— r 



Q„ = ^2B(n,r)-n r 

r=0 



and 



"' n-r 



F n = Y,^-rB(n,r) 2 u r+1 v n - r . 



r=0 r + 1 

To prove this one assumes a general form for the symmetric polynomial Q and adjust the 
coefficients so that the t-derivative is a total x-derivative. These functions also satisfy the 
equations 

t — ,a ^) 

J~ n,v ^=iLn,u j 

from which follows the conservation law (|J). The hierarchy of conservation laws may be 
combined into a single generating function 

so that the coefficient of A n in its formal power series expansion is the conservation law 
(|4]). To show this it is first necessary to find a recursion relation amongst the above charge 
densities. 

From the fact that Q n and T n are homogeneous polynomials in u and v of degree n 
and (n + 1) respectively one obtains the following relations 



nQ n = uQn >u + vQ n>v , , , 

(n + \)T n = uT n .u + vT n>v 

and 

F n +i,u = (n + l)(^n + vQ n ) , , v 

^n+i,, = (n + l)(^„ + «S n ). 1 ' 

Using these, together with (13), one easily obtains the following recursion relations: 



Q n+1 = ITn + [U + V) Q n , 

(n + 2)^ n+1 = (n + l)[(u + v)F n + 2uvQ ri 

which results in the second order recurrence relation for the charge densities: 

(n + 2)Q n+2 - (2n + 3)(w + v)Q n+1 + (n + l)(u - vfQ n = . 

This bears a strong resemblance to the recursion relation for Legendre polynomials. In 
fact, on defining 



Or. 



(u - v) n P n 



u + v 



u — V 



one obtains the Legendre recursion relation: 

{n + l)P n+ x{z) - (2n + l)zP n (z) + nP n - X (z) = 

and the first two conservations laws give P\(z) = 1 and P?,(z) = (3z 2 — l)/2 , the first two 
Legendre polynomials. Thus the charge densities may be expressed in terms of Legendre 
polynomials. Using the well-known generating function for these polynomials results in a 
generating function for the conservation laws: 



d_ 

dt 



u — v) 2 \ 2 — 2(u + v)X + lj 



d_ 

dx 



X(u + v) 



1 



u-v) 2 \ 2 -2(u + v)\ + l 



The form of this results suggest that one should define a new set of variables 



a u + v 
cos a = . 

u — v 

In this systems the equation {uQ u ) u = (vQ v ) v for the charge density becomes the axially 
symmetric Laplace equation in three dimensions. In terms of the inverse transformation 



u(x,t) = — r(x,t) cos 



2 0(x,t) 



2 



v[x, t) = r[x,t) sm — - — 



the original system (^) becomes 



n 

z t 



<r 2 /2)z x , 
(1 - z 2 )/2 r x 



(8) 



where z = cos 9 . The form of this suggest various generalizations, which will be investi- 
gated in section 5 . 

It is clear that any solution to axially symmetric Laplace equation, besides the class 
of solutions solutions r n Q n (cos6) will yield conservation laws. In particular one also has 
three other basic solution, namely 



-71—1 



PJ COS! 



-71—1 



Q n (cos6>) 



r n Q n (cosi 



where Q n {z) is the Legendre function of the second kind. If one wishes for conservation 
laws which remain regular as r — > (i.e in the limit in which the dispersionless Toda sys- 
tem becomes the dispersionless KdV equation) then this excludes the first two. However 
the third remains and gives a family of non-polynomial conservation laws for (|3|), the first 
few being 



1 1 t u \ 

2 iog y 



1 u 

-Qilog(-) - (u-v) 

2 v 



2S 2 log(-)--(« 2 



J t 



J t 



-(v — u) 



1 u 1 

-J zr 1 log(-) + ~(v-u)(v + u) 



1 u 1 

-J= 2 log(-) + -(v- u)(v 2 + lOuv + v 2 ) 

2 v o 



(the logarithmic terms in the (n + l) th such conservation law are of the form \ Q n \og(u/v) 
and |jF n log(M/f) ). Note that these are all antisymmetric under the interchange u *-> v 
and so vanish in the limit r — > . 



3. Generalised symmetries and commuting flows 

One fundamental property of an integrable system is the existence of an infinite family 
of commuting flows. For the system (||) it is straightforward to show that any such flow 
must be given by the system 

u t , = J r u x + uQv x , 
v t > = Tv x + v Qu x 

where T and Q are any function of u and v (not necessarily symmetric) which satisfy the 
relations @. Thus for each of the polynomial conservation laws derived in section 2 one 
has an associated flow labelled by the integer n (= degree of Q n ), 

Ut n = F n u x + u Q n v x , , v 

Vtn = FnVx+vQnUx- 



It is convenient to let to = t , = and Qo = 1 , so the original system corresponds to 
the to-flow. It follows from the general theory of generalised symmetries that the function 
J-u x + uQv x and Fvx + vQitx are characteristics for the original system, and that all these 
flows commute [HI . 



The charges are conserved with respect to the t -How, and one might expect them also 
to be conserved with respect to the t n -flow. This is indeed the case: 

Lemma 

9Qm =Fn 9Q m | Qn dF m 
dt n n dx n dx 
Furthermore, the right hand side is a derivative; 

dQ m _ dA(m,n) 
dt n dx 

where 

(m,n) = — — . (11) 

[n + l)[m + n + 1) 

Proof 

The first form of the equation (TiO| ) follows straightforwardly from (|9]). The second is 
most easily seen by the introduction of a potential. 

From (0) it is easy to change variables to obtain the Cauchy-Riemann like equations 

dF n dQ n 



dF m dQ m 
dQm 9F m 

Introducing a potential 



dA(m, n) dA(m,n) 



dQm ' dF m 
the equation (flO|) simply becomes 

dQm _ dA(m,n) dQ m | dA(m,n) dF m _ dA(m,n) 
dt n dQm dx dF m dx dx 

A more detailed analysis shows that the potential A(m,n) is given by ([11]). This makes 
use of equations (f|||) and (0). In the case m = 1, T\ — S, Qi — P the potential is simply 

given by A(l, n) = — 

K J n + 1 



4. The Hamiltonian structure of the hierarchy 



A system (||]) of hydro dynamic type is said to be Hamiltonian if there exists a Hamiltonian 
H = J dx h(u) and a Hamiltonian operator 

A^=g^{u) d - + U\{u)u k x 
which defines a skew-symmetric Poisson bracket on functionals 

^ ' ^ I 5u i (x) " 5ui{x) 
which satisfies the Jacobi identity and which generates the system 

ul = {u\x),H}. 

Dubrovin and Novikov || proved necessary and sufficient conditions for A 1 ^ to be a Hamil- 
tonian operator in the case when g l i is not degenerate. These are: 

a) g = {g 1 -')~ 1 defines a Riemannian metric, 

b) \) \ = —g ls Y 3 sk , where T J sk is the Christoffel symbol generated by g , 

c) the Riemann curvature tensor of g vanishes. 
The system (||]) may then be written as 

u\ = (g is V s Vjh)ui (13) 

where V is the covariant derivative generated by g . Thus to find the Hamiltonian structure 
for the hierarchy (^) one needs to find the metric and the Hamiltonians. 

One approach is to first diagonalise the system, there being the simple formulae for 
the metric coefficients for diagonal systems JJJ: 

djVj L, 

where V-(\i) = Vj(u)Sj (no sum). This diagonalisation is achieved by the transformation: 

U = v / w + y/v , 
V = \fu— \/v . 

While this enables one to find the metric, the formulae are actually more elegant in the 
{u,v} coordinates. 

Theorem 

The Hamiltonian structure is given by the zero-curvature metric 

du dv 
g = 2 

u v 

and the Hamiltonians are given by 



H n = J dxh n (u) 

where h n — (n + l) _2 Q n+ i . 
Proof 

The metric is clearly non-degenerate and flat, and the only non-zero Christoffel sym- 
bols are r 1 11 = — u^ 1 and T 2 22 = — v^ 1 (here u 1 = u and u 2 = v ). Hence, on expanding 
(P~3| ) one obtains the following equations for h n 



d 2 h n 

■Fn = UV TTTT ' 
ouov 

Q = 9 \ v dhn \ 

dv dv ' 

d dh n 
u- 



du du 

With h n — (n+ l)~ 2 Q„,+i these reduce to identities derived in section 2 and 3. Hence the 
result. 



This show the elegance of the {u, v} coordinate system: in terms of {P, S} or {r, 8} 
the metric is considerably more complicated. 

Proposition 

The Hamiltonians H n are time independent and in involution with respect to the 
Poisson bracket given by g . 

Proof 

The result follows from the formulae already derived. 



dH n , , 1 ._ 2 f , dQ n+ i 



(n + l) 2 J dx 



dA(n + 1, m) 



= o, 

under suitable boundary conditions. The second results follows from the use of the same 
formula: 



{H m , H n } = J 



dx dKn A ,dh n 



du 1 



(n + 1) 



J dxT n 



dQm+i T 9J-, m+ \ 

t" J~ 'n 



(n + 1 



i-2 



dx 



dx dx 
<9A(m + 1, n) 

dx 



The factor (n + l) 2 may, without loss of generality, be absorbed into a redefinition of 
the times t n . One therefore obtains the following generating function for the Hamiltonians: 



H{\) = J dxh(X) 



where 

MA) 



\^J(u-v) 2 \ 2 -2(u + v)\+ 1 A 



5. Generalizations 

The appearance of Legrendre polynomials suggests that there might be similar systems 
with conservation laws given in terms of more general polynomials. The system 



St — Px ) 

Pt = p w s x 



may be transformed into the system 



Ut = (uv) a uv. 
v t = (uv) a vu_ 



where a = (u — l)/(2 — u) (the case uj = 2 having to be treated separately). By an 
analogous proceedure to that outlined above, one can show that this has conservation 
laws expressible in terms of Gegenbauer polynomials , so when a = these collapse 
to Legendre polynomials. The Hamiltonian structure is also similar; this is given by the 
metric 

du dv 
g= u a+1 v a+1 ' 

and the Hamiltonians are directly related to the charges, as before. 

More generally still, one can look for polynomial conservation laws for the system 



r t = r a q(z)z x , 
z t = r b p(z)r x , 

this being a direct generalisation of equation (|8]). Looking for charge densities and fluxes 
of the form T = r c w(z) , Q = r d y(z) results in the constraints a — b = 2 and c = b + d—1 , 
and the self-adjoint differential equation 

d -( p { z ) d f)-d{d + a-l)y = Q (14) 



for y(z) . Assuming that this has regular points at , 1 and oo results in a hypergeometric 
equation and constraints on the otherwise free functions p(z) and q(z) . To get polynomials 
solutions requires d to be a non-negative integer, and so finally the charge densities are 

Qn = r n 2 Fi(-n, n + a - 1, c; z) , 
where 2 -Pi is the hypergeometric series. In terms of the u, v variables one obtains: 

v t = u a ~ c ~ 1 v c u x . 

There is, however, considerable degeneracy in this system, with apparently different sys- 
tems being connected via a change in variable. Other systems may be obtained, for 
example, by assuming that flUp is a confluent hypergeometric function, and this contains 
a family of charges given by Laguerre polynomials. 
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Appendix 

The results in this paper have been obtained by direct calculation; no use has been made 
of any Lax equation. In this appendix the Lax equation will be studied, and the solution 
to the associated linear problem derived. The system ([!]) is a special case of the more 
general equation 



d t d t , log P = d 2 x P 



known as the Boyer-Finley or SU(oo) Toda equation. This is itself a reduction of the 
equations for an anti-self-dual Ricci-flat metric. Physically the system (|IJ) describes anti- 
self-dual Ricci-flat metrics with two commuting Killing vectors, and so any Lax pair will 
fall into the class studies by Tod and Ward [IT]. One possible Lax pair for the Boyer- 



Finley equation (there are other, equivalent, formulations ||12|| ) was given by Ward |TJ 
and imposing the symmetry dt = df gives a Lax pair for (UJ). 

Explicitly, a Lax pair for the dispersionless Toda system (|TJ) is given by 



£0 = dt — Lf — £L e + , 
& = £(d t + L f )+L e -, 

where [Lf, L g ] = L {/i9} with {/, g} = f y g x -f x g y and e ± = Pt(x, t) exp(±y) , / = ±S(x, t) . 
Then the integrability condition [£ , C±] = for this otherwise overdetermined linear 
system 



£ v = o, 
Ci® = o, 

gives, on equating the coefficients of the various powers of £ , the spectral parameter, to 
zero, the system P t = PS X , S t — P x . 

The solution \P to the linear system may be found in the form of a power series 

oo 

* = J2^n(x,t) exp(ra/)f\ 

n=0 

The if} n are not unique - one can multiple by an arbitrary holomorphic function of £ and 
this will change all the coefficients, but not the integrability conditions, this corresponding 
to a gauge freedom. This yields an infinite number of recursion relations 



d x 9 n+ i — P n d t (P n 9 n ) , 

d t 9 n+1 = P n+1 d x (p- n 9 n ) , 

where 

B n (x,t) = 4> n P 2 ■ 

together with the initial conditions d x 8 = d t 9\ = . The first few solutions to these 
equations are 



9 = k , 

6 1 = h, 

6 2 = -k 1 S + k 2 , 

3 = h(S 2 - P) - 2k 2 S + k 3 . 

Here the ki are constants. Using the above-mentioned gauge freedom one can set ko = 
and ki — for % > 2 , so that, for n > 1 , 9 n is an (n — l) th polynominal in u and v . After 
a few more iterations one finds the relation 

On+i + S6 n + PQ n _ x = , n > 2 

appears to work, and this can be proved by induction (again using the the gauge freedom 
to set various constants to zero). Solving this relation gives 

/ „,n „,n\ 

n+lj 



e n = (-i) n+L h 



u — V 



and hence ip n . Summing the power series is straightforward: 

fci 1 



£ 1 exp(-y) 



£exp(+y) 



where £ = £ exp y . In all this y plays an auxiliary role. This may be removed in favour of 
£ = £ exp (+y), but this will introduce terms involving in the Lax pair. 
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